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Even though the rare-earth tritellurides are tetragonal materials with a quasi two dimensional 
(2D) band structure, they have a "hidden" ID character. The resultant near-perfect nesting of the 
Fermi surface leads to the formation of a charge density wave (CDW) state. We show that for this 
band structure, there are two possible ordered phases: A bidirectional "checkerboard" state would 
occur if the CDW transition temperature were sufficiently low, whereas a unidirectional "striped" 
state, consistent with what is observed in experiment, is favored when the transition temperature is 
higher. This result may also give some insight into why, in more strongly correlated systems, such 
as the cuprates and nickelates, the observed charge ordered states are generally stripes as opposed 
to checkerboards. 



I. INTRODUCTION 

Many strongly correlated (layered) materials with 
quasi two dimensional (2D) structure show fluctuating 
or static stripe order (unidirectional density wave states) 
over a substantial range of temperaturesiSi^ In par- 
ticular, there has been considerable interest in study- 
ing whether or not stripes are inextricably connected to 
the high temperature superconductivity that occurs in 
the cupratesi^ However, characterizing the stripe order 
in these materials is difficult due to the complex and 
strongly correlated nature of these materials, and the 
presence of significant amounts of quenched disorder^ 

From the symmetry view point, there is no differ- 
ence between a unidirectional charge density wave state 
(CDW) in a weakly interacting quasi-2D system and 
a stripe-ordered state of a strongly correlated system. 
CDW states can occur in the weak coupling limit only 
if there are sufficiently well nested portions of the Fermi 
surface. Since this is non-generic in more than ID, one 
would like to identify what is special about those higher 
dimensional materials in which nested Fermi surfaces ap- 
pear. Moreover, for a layered quasi-2D material with 
tetragonal (C4) symmetry, the CDW ground state can 
either be bidirectional (checkerboards) maintaining the 
point group symmetry of the lattice, or unidirectional 
(stripes) with a reduced symmetry, and again we would 
like to understand what physics governs this choice. 

The rare-earth tritclluride series i?Te3 (R=Y, La-Tm) 
is particularly well suited for a detailed study of these is- 
sues. i?Te3 consists of square Te planes alternating with 
weakly coupled i?Te slabs. (See Fig. 2 of Ref.0 for exam- 
ple.) The electronically active valence band derives from 
the 5p orbitals of the planar Te atoms, and is thus ex- 
pected to be relatively weakly correlated. The existence 
of a unidirectional CDW was first detected by trans- 
mission electron microscopy (TEM)Ai Angle resolved 
photoemission spectroscopy (ARPES)^iS measurements 
have shown that the CDW ordering wavevector nests 
large portions of the Fermi surface, which are in turn 
gapped, thus strongly indicating that the CDW is a con- 
sequence of Fermi surface nesting. More recent X-ray 
scattering and STM measurements have confirmed the 



existence of the CDW and its unidirectional character 
in great detail>i2*ii Detailed studies of the Fermi surface 
topology using ARPES^ and positron annihilation^ sup- 
port the notion that the CDW is the result of a nesting- 
driven Fermi surface instability. However, the driving 
force for the strong breaking of the point group sym- 
metry (C4 — > C2) produced by the unidirectional CDW 
formation in -RTC3 has not been clear previously. 

We shall show in the present paper that the nest- 
ing of the Fermi surface of -RTC3 reflects a "hidden" ID 
character of the electronic band structure which derives 
from the highly anisotropic hopping amplitude of the 
Te p x and p y orbitals^ Consequently, CDW order oc- 
curs for a dimensionless effective interaction, A, in ex- 
cess of an extremely small critical value, A c « 1. (In 
ID, the fact that A c = is the famous Peierls instabil- 
ity.) Moreover, we shall show that there are two pos- 
sible patterns of CDW order that can take best advan- 
tage of the nesting: checkerboard order that occurs for 
A slightly larger than A c , and stripe order, which is also 
rotated relative to the checkerboard order, which occurs 
for A > A ~ A c [l + C(A C )]; in more physical terms, what 
this means is that when the CDW transition tempera- 
ture, T c , or the CDW gap, Aq, is sufficiently large, stripe 
order is favored. (The resulting phase diagram is shown 
in Fig. 4.) 



II. BAND STRUCTURE CONSIDERATIONS 

It is well established from transport measurementsi^i^ 
that the coupling between layers in i?Te3 is small. Com- 
parisons between different rare earth compounds con- 
firm the minor role of the rare earth in the electronic 
structureAi^ Hence the physical properties of i?Te3 are 
dominantly determined by the Te planes common to all 
i?Te3. It is therefore natural to consider a simple model 
of the electronic structure of a single Tc plane. 

It is known, from first principle band structure 
calculations^ that the energy of the hp z is shifted by 
crystal field effects so that it lies well below the Fermi 
energy. Consequently, the relevant bands close to the 
Fermi surface are the 5p x and 5p y orbitals, which arc 




FIG. 1: The p x (hatched) and p y (solid) orbitals in the sin- 
gle Te square lattice. Due to the highly anisotropic profile of 
the p orbitals, the hopping amplitude, t\\ , along the extended 
direction of the given p orbital is much larger than the hop- 
ping amplitude t± along the direction perpendicular to the 
extended direction, if is a second-neighbor hopping, which 
is the shortest-range interaction which mixes the p x and p y 
bands, and which we neglect for present purposes. 
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FIG. 2: Fermi surfaces (a) for t± = and (b) for t±_ small, 
but non-zero. For t±_ — 0, the wavevector Q2 = (2kF,2kp) 
nests both the p^-FS and the p y -FS perfectly. For non-zero 
tj_, the wavevector Qi = (2kF,ir) does a much better job of 
nesting the Pr-FS but a poor job on the p y -FS. On the other 
hand, Q2 does a moderately good job of nesting the entire 
FS. 



approximately 5/8 filled. The crystallographic unit cell 
contains two Te atoms per plane due to a slight in- 
equivalence between two sublattices produced by the 
i?Te slab orientation and the crystal structure is very 
weakly orthorhombic^ for simplicity, we will ignore 
these subtleties and consider an idealized Te square lat- 
tice. The band structure can be well approximated in the 
tight-binding approximation; the resulting Hamiltonian 
is shown schematically in Fig. ^ 

Following Ref. 0, we neglect all interactions of longer- 
range than the nearest neighbor hopping, with the con- 
sequence that there is no hybridization between the p x 
and p y bands. (Small further neighbor interactions, such 
as the second neighbor hopping t' in Fig. ^ have a large 
effect on the Fermi surface topology where the two bands 
intersect, but can readily be shown to have little effect on 
the results obtained below.) Working in units where the 
lattice constant of the square lattice a = l, the dispersion 
for the p x band and p y band can be readily derived: 



£k,p x = — 2£|| cos(k x ) + 2t± cos(ky) 
\(k x ) - 2t|| cos(ky), 



£k,p B = 2i j_ cos 



(2.1) 



where tu (tj_) represents the hopping amplitude parallel 
(perpendicular) to the extended direction of the given 
p orbital. Due to the highly anisotropic profile of the 
p orbital electron wave function, the hopping amplitude 
tu , along the extended direction of the given p orbital, is 
much larger than the hopping amplitude t± perpendic- 
ular to the extended direction. Indeed these hopping 
amplitudes have been estimated^ to be tu « 2.0eV, 
t± « 0.37ey, and if w 0.16cV for i?Te 3 . Thus, it is 
reasonable to set t' = and to treat tj_/tu « 0.18 as a 
small parameter. 

The small magnitude of the ratio t±/tn implies a se- 
cret quasi-lD character of the band structure. For t± = 0, 
the system would be equivalent to an array of ID wires. 



Nonetheless, the system would maintain overall C4 sym- 
metry since these p x and p y "wires" are perpendicular to 
each other. Thus, even in this limit, the system would 
not display 2D anisotropy in any transport measurement. 
However, the resulting band structure would consist of 
two parallel ID FS's as shown in Fig. Ela). With kp 
defined by the implicit relation p = — 2tu cos kp, where 
(i is the chemical potential, any wavevector (±2kp, k y ) 
would then perfectly nest the p^-FS for arbitrary k y and 
wavevectors (k x ,±2kp) would perfectly nest the p y -FS 
for any k x . In particular, the wavevectors (±2fcp, ±2kp) 
perfectly nest both the p x and p y FS's. 

A small, but nonzero t±_ introduces nonzero but small 
curvature to the FS (see Fig- I2tb)); the p^-FS is deter- 
mined by 
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(2.2) 



where kp is defined above. Here +(— ) represents the 
right (left) portion of the p^-FS respectively. Obviously, 
the Py-FS can be obtained from the p^-FS by 90° rota- 
tion. 

As a consequence of the finite curvature, the wavevec- 
tor (±2kp 1 k y ) nearly perfectly nests the p^-FS only for 
ky = 7T, and equivalently the p. y -FS is best nested by 
(it, ±2kp). However, away from the half filling, (±2/cp, 7r) 
does a poor job in nesting the p y -FS and conversely for 
(T7,±2k F ). Note that (2fcp,7r), (-2kp,ir), (tt, 2k F ) and 
(■7T, —2k F ) are related by the C4 symmetry of the host 
lattice. On the other hand, the wavevector (2/cf,2/cf) 
(and its C4 symmetry related vectors), which perfectly 
nests the full FS for t± = 0, does a reasonable job in 
nesting both the p^-FS and the p y -FS. Thus, for small 
there arc two independent candidates (not related 
to each other via C4 symmetry) for the CDW ordering 
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wavevectors: 



Qi = (2fc F;7 r) and Q 2 = (2k F , 2k F ). (2.3) 

To sec which ordering vector is preferred, we now com- 
pute the corresponding charge density (Lindhard) sus- 
ceptibility, x(<l',T). For a given wavevector q at a given 
temperature T, 



x(q;T) = -2 



v 



dk f(& v ,T)-f(& +m ,T) 



(2-Kf 



6c i) - 6 



k+q?) 



(2.4) 

where the integral is over the first Brillouin zone, rj = 
p x ,Py is the band index, and £kr; = £kr; — P- The factor 
of 2 is from two spin polarizations of the electrons and 
K^krpT) denotes the Fermi-Dirac distribution function. 

For small t±/t\\, the susceptibilities at these two vec- 
tors have the following approximate analytical forms, 
which are valid for temperatures T <C t± (in units 
k B = 1): 



x(Qv,T) n p(E F )log 

+ p(E F )k 
X (Q 2 ;T)«2p(£ F )log 



a t\\ 



T+ ai t\m 

a t\\ 



T + a>2t\\ 
aot\\ 



T + a 3 t± 



, (2.5) 
(2.6) 



where p(E F ) = [27rtii sinfc F ] _1 is the density of states per 
spin per band and «o^|| is the ultraviolet cutoff, which is 
of order of the bandwidth (a ps 2). The parameters oti 
arc defined by 



a\ ~ | cos k F 1/(4 sin 2 k F ), 
a 2 w 2 1 cos k F | , 
a 3 k, | cos k F | . 



(2.7) 



In i?Te3, the bands arc approximately 5/8 filled, thus 
k F w 5n/8. In the small but non-zero t±/t\\ limit, 

x(Qi;0)» x(Q 2 ;0) + p(i? F )log[2sin 2 (A :F )]. (2.8) 

Therefore, if the CDW transition occurs at low enough 
temperatures the wavevector Qi will be favored so long 
as 7r/4 < kp < 3ir/4. However, at high enough tem- 
peratures, T > T ~ t\/t\\, x(Q2',T) is greater than 
x(Qi;T). Therefore, we expect there to be a finite tem- 
perature To such that if the CDW transition temperature 
T c <Tq, the ordering vector is q w Qi , while for T c >Tq, 

q ~ Q2 

For a more quantitative analysis, we have numeri- 
cally evaluated the susceptibility x(q; T) of Eq. I|2.4|l . for 
wavevectors q over the whole Brillouin zone. The calcu- 
lated susceptibility has a maximum at a wavevector q max 
close to Qi or Q 2 , depending on whether the tempera- 
ture, T < To or T > To, in agreement with our analysis 
of the FS. The precise location of (\ max shifts slightly as 
a function of T, but so long as T C tii, the maximum 
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FIG. 3: The Linhard susceptibility x(cii T) as a function of T 
for q = Qi = (2k F , n) and q = Q2 = (2k F , 2k F ). 



always lies close to Q ^ , so this deviation does not signif- 
icantly affect the resulting picture of the phase diagram. 
We will henceforth neglect this small effect, and focus 
our discussion on the wavevectors Qj. Fig. [3] shows the 
calculated x(<i\T) for q=Q^. 

This result can be understood as follows: At low tem- 
peratures, the states arbitrarily close to the nearly per- 
fectly nested portions of the Fermi surface dominate the 
susceptibility favoring Qj. However, for T > T , the cur- 
vature of the FS has negligible effect, so Q2 is preferred. 
Since Qi nests only the p x portion of the FS, there is lit- 
tle competition between the ordering tendency at Qi and 
at the symmetry related vector, Qi = (tt, 2k F ). Thus, 
as we shall show below, wherever CDW ordering with 
wavevector Qi occurs, simultaneous ordering occurs at 
Qi = (ir, 2k F ), making this a state with "checkerboard" 
order. However, since ordering at wavevector Q 2 opens 
a gap on the entire FS, it prevents ordering at the sym- 
metry related vector, Q 2 = (—2k F ,2k F ); when Q 2 is the 
preferred ordering vector, the resulting state has "stripe" 
order. 



III. MEAN-FIELD PHASE DIAGRAM 

To obtain an explicit phase diagram, we add an 
electron-phonon coupling and solve the resulting model 
in mean-field approximation. Thus, we investigate the 
following mean-field Hamiltonian 



H 



MF 



k,r) q 



p(Ef 
2\„ 



|Aq| 2 



+ E ^k+q^k,, (3.1) 
k,q,») 

where the sum over q runs over the possible ordering vec- 
tors, q = ±Qi, ±Qi, ±Q 2 , and±Q 2 , and all harmonics, 
such as ±2Qj and ±Qj ± Qj. The order parameter A q 
for the CDW with the ordering vector q is related to 
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the lattice distortion (x q ) through A q = a q (:z; q ), where 
a q is the electron-phonon coupling, and the dimension- 
less effective interaction A q = [a^p(Ep)]/[Muj'^\, where 
w q is the phonon frequency and M is the ion mass. The 
self-consistent value of the gap is obtained by minimizing 
the free energy computed from Hmf with respect to A q . 
Note that A_ q = A*. We use the approximation A q = A 
in the rest of this paper. 

At high temperatures, clearly A q = for all q. As the 
temperature is lowered, for sufficiently large A, solutions 
with non-zero values of A q appear. In general, when 
A q is non-zero at some wavevector, q, it is non-zero, 
although possibly much smaller, at all harmonics, as well. 

For commensurate order, the mean- field equations can, 
in principle, be evaluated numerically exactly, but for 
high order commcnsurability this is quite a complicated 
problem, as there are many harmonics, and hence many 
possible patterns of order. For incommensurate order, 
even this is not possible. 

Fortunately, in the present case, for the most part the 
transitions are continuous, or at worst weakly first or- 
der. Therefore, the phase boundaries can be identified 
by expanding the mean-field (Landau) free energy to low 
order in powers of A q . The coefficients in this expansion 
are expressed as convolutions of free Green functions, as 
derived explicitly up to fourth order in the Appendix. 
To this order, only the fundamentals and their second 
harmonics enter the expansion; if q and q are two (sym- 
metry related) fundamental ordering vectors, then the 
higher harmonics in powers of the fundamental, satisfy 
A„ q+mq ~ [AqAq 1 ]. Indeed, in this context, we can treat 
the harmonics as slaves to the fundamentals and obtain 
a Landau free energy expressed exclusively in terms of 
A q and A q , as is done explicitly in the Appendix. The 
result is a free energy density 



2 

^ = ^[r J (|A Qj .| 2 + |A Q J 2 ) (3.2) 




where 

r j =p(E F )/X- X (Q j ;T), (3.3) 

and . . . signifies both higher order terms in powers of 
A, as well as unimportant biquadratic terms, such as 
| Aq 1 | 2 | Aq 2 | 2 , which come into play only if ordering oc- 
curs both at Qi = (2kp,ir) and Q2 = (2fcp, 2fcp). The 
coefficients Uj and are evaluated in the Appendix. 

The mean-field phase diagram in Fig. 01 was con- 
structed from this expression as follows: The phase 
boundaries were determined as the points at which ei- 
ther ri or T2 vanishes, so both t\ and T2 > in the high 
temperature phase. Throughout the ordered region of the 
phase diagram, we always find that 71 < and 72 > 0, so 
wherever Qi ordering is favored, simultaneous ordering 
at Qi and Qi, i.e., checkerboard order, occurs, while in 
the regions where Q2 ordering is favored, either Aq 2 or 
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FIG. 4: The mean-field phase diagram. For A > Ao, the system 
orders to a stripe phase below the transition temperature. 
The associated phase transition is second order for A > A tc 
and first order for Ao < A < A tc - For A c < A < Ao, the system 
orders to a checkerboard pattern. For A < A c , there is no 
CDW instability down to zero temperature. Explicit numbers 
refer to the representative values of t\\ — 2.0eV, tj_ — 0.37eV 
and kF ~ 57r/8. Tmf signifies the mean field T c ; the actual 
T c will be smaller due to fluctuation effects. 

Aq 2 is non-zero, but not both, i.e., striped order occurs. 
The transition to the stripe phase is continuous so long 
as «2 > 0, which is the case at large coupling. The tri- 
critical point occurs where ui vanishes. Strictly speaking, 
where the transition becomes first order (the dashed lines 
in the figure), the Landau expansion is inadequate; the 
dashed lines should therefore be interpreted as schematic 
representations of the true phase boundaries. In evaluat- 
ing the coefficients in the Landau expansion, we use the 
band structure with the representative values given in the 
caption of Fig. but the phase diagram is qualitatively 
similar for any t± <C tn and 7r/4 < kp < 3ir/4. 

The main new physics apparent in this phase diagram 
arises from the sign of jj. The fact that 72 is robustly 
positive follows from the fact that Q2 does a moderately 
good job of nesting the entire FS; in the presence of a 
non-zero Aq 2 , most of the FS is gapped, so there are 
no remaining portions of the FS to be nested by Q2- 
Conversely, the fact that Qi fails so badly to nest the p y 
portions of the FS, implies that there is no substantial 
interference between Qi and Qi ordering; this leads to 
the expectation of a near vanishing biquadratic coupling, 
or in other words, 71 « —u\/2 < 0, an expectation that 
we find is approximately satisfied. (At T = 0, 71 « —25 
and ui « 47.) 

Using the Landau free energy with coefficients given in 
the Appendix, for the band structure with the representa- 
tive values, we can estimate the various critical values of 
the of the coupling constant: the minimum value for the 
existence of a CDW ground state is A c « 0.093; the criti- 
cal value that separates the regime of checkerboard order 
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from stripe order is Ao ~ 0.103; the tricritical coupling 
is A tc « 0.105; the temperature at the tricritical point 
is T tc w 0.098 cV. The maximum checkerboard ordering 
temperature, To, which occurs on the edge of the first or- 
der transition to the striped phase, is the one qualitative 
feature of the phase diagram that depends on the higher 
order terms in the Landau expansion that we have not 
computed. The estimate for To given in the figure comes 
from the assumption that the first order phase boundary 
lies close to the point at which r 2 vanishes. 



IV. FLUCTUATION EFFECTS 

Because of the hidden ID character of the ordering we 
have been exploring, we expect fluctuation effects will 
have a large quantitative effect on the mean-field phase 
diagram. Indeed, it is common in quasi ID CDW sys- 
tems for the ordering temperature, T c to be significantly 
suppressed below its mean-field value. This is reflected in 
larger values (than the prediction of mean-field theory) 
of the ratio, Ao/T c , of the zero temperature gap to the 
actual T c p±L±£ 

We can estimate the extent to which fluctuations sup- 
press T c in two different limits: If the anisotropy is 
not too large, i.e. if t\\ 3> t± ^ Ao, the suppres- 
sion is fractionally small, and can be estimated, as in 
the theory of superconductivity, by the Brout criterion, 
Tmf ~T c ~ T MF Almti_. Conversely, if t\\ > A > t±, 
while Ao may still be crudely determined by mean-field 
considerations, T c <C Tmf an d is determined by entirely 
different physics. In this limit, the intra-chain CDW fluc- 
tuations can be treated using the theory of the one di- 
mensional electron gas, and t± can be included in the 
context of inter-chain mean-field theory, with the result 
that T c ~ t±(tj_/t^) a where a is an interaction depen- 
dent constant Whenever there is a large fluctua- 
tional suppression of T c , local CDW correlations are ex- 
pected to survive in a broad range of temperatures above 
T c ; roughly, the local CDW correlations develops in the 
temperature range T c < T < Tmf, where Tmf is the 
mean-field transition temperature. 

For i?Tc3, Ao ~ 260 — 400mcVj£ corresponding to a 
mean transition temperature Tmf ~ 1500 — 2000K and 
A ~ 0.1. Since t± sa 0.38cV ~ Ao, a reliable estimate of 
T c is not possible for these materials, although a large 
supression relative to the mean-field value can be ex- 
pected. An expression due to Barisio^i based on physi- 
cally plausible but hard to justify approximations yields 
T c ~ {t±/t\\)TMF ~ 300K. Indeed, very recently, it has 
been found^ that he CDW phase transition in i?Te3 oc- 
curs (depending on R) in the range T c = 260 — 400/, K, 
and substantial CDW correlations persist well above T c . 

One consequence of remaining local CDW correlations 
above T c will be the presence of peaks in the structure 
factor at positions corresponding to the Bragg vectors 
of the ordered state, but with width inversely propor- 
tional to the thermal correlation length^ However in 



the case of a stripe phase, more dramatic effects can 
be expected. If there is only a single transition; then 
the "stripe liquid" phase above T c does not break any 
of the lattice symmetries, and consequently equivalent 
peaks in the structure factor should appear both at the 
stripe ordering vector, and at the conjugate wavevector 
(rotated by 90°). Such behavior has already been seen, 
albeit only in the magnetic scattering, in the stripe liquid 
phase of La2- x Sr x Ni04i24 On the other hand, a two stage 
transition is also possible with an intermediate, "stripe 
ncmatic" phased in which stripe fluctuations are suffi- 
ciently violent to restore translation, but not the full C4 
symmetry of the host crystal; in such a phase, the peaks 
at the stripe ordering wave vector should be stronger 
(possibly, much stronger) than those at the conjugate 
wavevector. In this case, there must be a second transi- 
tion at still higher temperatures to an isotropic state. 

If future experiments can confirm the predictions above 
regarding the fluctuation effect, they will provide an im- 
portant laboratory for exploring the physics of a stripe 
liquid, with possible broader implications to many stripy 
materials. 



V. CONCLUDING REMARKS 

The present mean-field theory gives results that are 
broadly consistent with experimental observations in 
i?Te3. From the value of the CDW gap measured by 
ARPES)^ one can estimate the mean- field transition 
temperature to be Tmf ~ 2000 K, which is an order 
of magnitude greater than the theoretical value of To 
and also greater than T tc . Hence, the theory predicts 
that i?Te3 should have a unidirectional CDW ordered 
phase (stripes) with ordering vector approximately equal 
to cither Q 2 = (2k F ,2k F ) or Q 2 = (2k F ,-2k F )& This 
is consistent with the experiments^*^ which find stripe 
order with Qcdw ~ (0.717T, 0.7l7r) in the present non- 
doubled unit cell notation^ With the simplified band- 
structure used in the present analysis, 2k F w 37r/4 (mod 
27r), However, as mentioned above, the maximum of the 
susceptibility occurs at a slightly different weakly tem- 
perature dependent wave-vector; at T = 350K, the max- 
imum occurs at q mQX = (0.72-7T, 0.727r), which is almost 
identical to the experimental results at comparable tem- 
peratures. Moreover, the mean field theory predicts that 
the phase transition into the stripe ordered phase from 
higher temperature phase is a second order phase tran- 
sition. This is also consistent with recent experimental 
results^ 

Even though we are dealing with a stoichiometric phase 
with an integer number of electrons per unit cell, if the 
transition is continuous, then (at least for commensura- 
bility N > 4) the CDW ordered state just below T c is 
generically incommensurate, in the sense that the order- 
ing vector is not locked to the lattice. This follows from 
the fact that the lowest order term which locks the CDW 
to the lattice is proportional to A N cos(N(f>) and so is ir- 
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relevant (for N > 4) in the renormalization group sensed 
at T c . (Here A is the order parameter and (f> = 2kpx$ 
determines the relative phase between the CDW and the 
underlying crystalline lattice.) In the present problem, 
two effects contribute to shift the ordering vector from 
its commensurate value. Firstly, even at T — 0, there 
is a shift of Uf = 37r/4[l + 0(t±/t\\) 2 ] (which we have 
not discussed explicitly) due to the 2D dispersion. More- 
over, Q2 depends weakly on the temperature and (as 
mentioned above) is noticeably different from its T = 
value at T = T c . However, it is possible that, under 
appropriate circumstances, there will be a second com- 
mensurability lock-in transition with critical temperature 
Tcom < T c . 

Another subtlety of the problem derives from the ex- 
istence of Te bilaycrs in -RTC3. If the bilayer split- 
ting, tbu <C Ao, it can be ignored, and a single tran- 
sition occurs as treated in the present work. However, 
if tit ^> tbu Ao, we should treat the ordering in the 
bilayer split bands separately. In this case, there should 
be two distinct ordering transitions, at distinct order- 
ing temperatures, with slightly different ordering vectors, 
kp,± = k F ± t b ii/v F . 

Recent results^ on RTe3 with different rare earth el- 
ements, R, show that in some cases there does, in fact, 
appear to be a second phase transition at temperatures 
below T c . For various reasons, we conjecture that where 
the second transition occurs, it is due to the bilayer split- 
ting, but whether it is this, or a commensurate lock-in, 
or some other transition remains unsettled. 

Charge density waves often break the point group sym- 
metry of the host lattice. This phenomenon is particu- 
larly common in strongly correlated materials including 
many transition metal oxides which exhibit stripe order. 
The strong interactions present in those (d-band) sys- 
tems make the physics of pattern selection more diffi- 
cult to study from a microscopic viewpoint. However, 
one may hope that extrapolating from the weak coupling 
limit (solved in the present paper) may give some insight 
into the deeper issues of pattern selection in highly cor- 
related materials. 

The underlying physics that is responsible for the ex- 
istence of a CDW state in the present class of reasonably 
weakly interacting quasi-2D systems is the existence 
of a hidden ID character of the band structure. In 
the present paper, we have explored only the grossest 
aspects of this structure - primarily those amenable to 
a mean- field analysis. At temperatures or frequencies 
in excess of the ordering temperature, where the band 
structure can be approximated as that of intersecting ID 
systems, it is probable that more interesting fluctuation 
effects, associated with the breakdown of Fermi liquid 
theory in ID, should be observable. Indeed, some evi- 
dence already exists^ from high energy spectroscopies 
of anomalous power-law behaviors reminiscent of the ID 
Luttinger liquid. We are currently exploring this aspect 
of the problem. 
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APPENDIX A: THE LINKED CLUSTER 
EXPANSION 
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FIG. 5: The second order Feynman diagram from which we 
obtain the Lindhard susceptibility, \Q- Higher order dia- 
grams differ only in the number of external legs. 

Here we sketch the details of the linked cluster expan- 
sion that was used to obtain the coefficients of the GL 
free energy. We treat the off-diagonal part of the mean 
field Hamiltonian Eq. (|3.1|) 



k,q,?7 



(Al) 



as a perturbation to obtain the Landau free energy: 
^[{A q }] = To + ^ 4^ |A q | 2 + f] t/ ; [{A q }], (A2) 



^ 2A 
q 



1=1 



where the set {A q } contains, in principle, all harmonics 
of the fundamental ordering vectors, although in practice 
we restrict ourselves to fourth order terms for the contin- 
uous phase transitions. The functionals, Ui, are defined 

as 



(-i) e+i 1 r p 







dn(T T V{n)---V(Tt)). 

(A3) 



The expansion in terms of U( contains only fully con- 
nected, 1-loop diagrams, where the finite-temperature 
Green functions for the electrons are 



<$(k,w n ) = 



1 



(A4) 



with uj n = (2n + l)ir/0, and the order parameters, Aq, 
act as external, single-particle potentials (see, for exam- 
ple, Fig. EJ). 

To fourth order, the Landau free energy in terms of 
order parameters associated with Qi, Q2 and their har- 
monics is 



G 



T = ^ + E {^(l A Q, I 2 + I A Q 3 1 2 ) + A 2 Qj I 2 + |A 2Qj | 2 ) + f, (|A Qj+Qj | 2 + |A Qi _ Q . | 2 ) + d,(|A Q3 1 4 + |A Qj | 4 ) 

3=1 

+ 9j \A Ql | 2 1 A Q2 1 2 + [6, (A 2 Q . A_ 2Qj + A| . A_ 2Q2 ) + Cj ( A Qj A Q . A 

Q Q A Q A Q A Q Q 

) + c.c.]},(A5) 

where the coefficients are given by the integrals 

r = p(E F )/X- X (Q 3 ;T) (A6) 
fj = p(E F )/X- x (2Qj;T) (A7) 
r, = p(E F )/X - x (Qj + Qj;T) (A8) 

h = /"e?(0)Sj(Q i )eJ(2Q i ) (A9) 



{^(0) Go(Qj)G%(Qj + Qj) + Sg(0) ^o(Qj) (Qj + Qj)} (Aio) 

<*J = /'{^(0)^(Q J ) 2 ^(2Q J ) + 5^(0) 2 ^(Q,) 2 } (All) 

{^(o) 2 ^(Q,)ao'(Q,) + ^(o) 2 ^(Q,)^(-Q.) + 6?(o)e?(Q 3 -) 2 e2(Qi + Qj) 

^'(0)^(Q,) 2 ^'(Q J - Qj) + 2Sg(0)Sg(Q 3 -)6S(Qi)6J(Qj + Qj)}, (A12) 



I 

and we have adopted the compactified notation: The effect of higher harmonics is to renormalize the co- 

efficients of the quartic terms, which become 



/ 



ej(q)=fij(k + q,iw n ) (A13) 



?7 za;^ 



As the susceptibility towards ordering is always largest 
for the fundamental vectors, we integrate out the higher 
harmonics, which allows us to write the free energy in a ^. = g _ \cj\ 2 /fj — uj/2. (A16) 

form that transparently displays the C4 symmetry of the 
problem 

^ = ^h(|A Qj | 2 + |A Q J 2 ) (A15) 

3 

+ ^(|A Q3 | 2 + |A Qj | 2 ) 2 + 7 ,|A Qj | 2 |A Q2 | 2 ]. 
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